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We study a novel S — 1/2 cluster chain Hamiltonian which has recently been proposed in the 
context of the charge ordered low-temperature phase of Q'-NaV205. We perform a detailed inves- 
tigation of this model within a large range of parameters using perturbation theory and Lanczos 
diagonalization. Using model-specific local conservation laws and parameter-dependent mappings 
to various effective low-energy Hamiltonians we uncover a rich phase diagram and several regimes 
of gapful spin-excitations. We find that the overall features of recent neutron scattering data on 
a'-NaV205 can be fitted within this model, however using a set of parameters which seems unlikely. 

PACS numbers; 75.10.Jm, 75.40.-s, 75.40.Mg, 75.50.Ee 
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I. INTRODUCTION 

Low-dimensional magnetic materials exhibit many in- 
teresting and sometimes puzzling properties. In partic- 
ular, the origin of the spin gap in the low-temperature 
phase of a'-NaV205 is under intense discussion since it 
was first reported to open with a phase transition at T = 
34K pi. It was confirmed soon [p| that this transition is 
accompanied with a lattice distortion. Charge ordering 
was also realized to play an important role |g|| , leading to 
various ordering patterns proposed M-M ■ Recently, based 
on new determinations of the low-temperature crystal 
structure JT) , two groups |g,g| have proposed a new type 
of charge ordering which is sketched in Fig. |l|. This struc- 
ture incorporates Vanadium ions in three types of valence 
states: V*"*", V^"'' and V*'^^. Entering the low-tempera- 
ture phase below T — 34K, every second of the quarter- 
filled two-leg ladders, which make up the system at high 
temperatures pO| ] is in a zig-zag ordered pattern proposed 
previously for all of the ladders [g] . The remaining lad- 
ders survive the transition. This new proposal does not 
completely rule out other scenarios e.g. with only V'*''' 
and V^^ since the charges were determined only indi- 
rectly through the crystal structure whose determination 
in itself appears to be delicate [Q. Nevertheless, the 
charge order proposed by ||,g| is appealing since it con- 
tains clusters of six Vanadium ions with a total of four 
unpaired electrons (indicated by bold lines in Fig. H) and 
therefore a spin gap arises naturally H. 

The nature of magnetic excitations above the gap in 
the low-temperature phase of a'-NaV205, i.e. the dis- 
persion of 5 = 1 excitations, has been determined by 
neutron scattering [n3-n4|: Along the a-direction, two 
weakly dispersive excitations with an approximate band- 
width of ImeV are observed. In contrast, there is a rather 
large dispersion along the fe-direction ||l^, with a band 
rising from around lOmeV to at least 40meV [H. This 
suggests that a one-dimensional model extending along 
the 6-axis should result in a reasonable first approxima- 
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FIG. 1. Possible charge order of Vanadium ions in 

a'-NaV205 (schematic). Lines denote the most important 
(super)exchange paths. The Vanadium atoms are located at 
intersections and corners; valence states are indicated by the 
numbers in the center of the figure. Grey shaded areas indi- 
cate the locafization regime of one unpaired electron. Decou- 
pling along the dashed lines leads to the quadrumer chain. 



tion of the magnetic properties. In the charge-ordered 
state of Fig. Ill, this can be implemented by a decou- 
pling along the dashed lines. To leading order in 1/U 
of the related Hubbard model, this leads to the spin-1/2 
quadrumer chain of Fig. [|l5| . 

The properties of the cluster chain model and its rela- 
tion to a'-NaV205 are subject to controversial discussion 
|l5| , |l6| . Moreover, not all regimes of parameters possi- 
bly relevant to a'-NaV2 05 have yet been considered in 
the literature. In particular, the case J' < remains 
to be investigated as some LDA-I-U calculations seem to 
indicate a negative J' ]1^. Apart from its possible rele- 
vance to a'-NaV205, the quadrumer model is closely re- 
lated to several other interesting spin-chain Hamiltonians 
|18|-pi|. Eg., CaV409 was discussed in the framework of 
the model in Fig. 0, however with J2 = [18|. The latter 
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FIG. 2. The quadrumer-chain Hamiltonian. S = 1/2 spins 
are located at the numbered intersection points and corners. 
Coupling constants are indicated by the style of the lines. 

model was also generalized [|l9| by adding vertical cou- 
plings in the quadrumers of Fig. 0. Local reflection sym- 
metry of these cluster spin models along the chain direc- 
tion allows for exact statements regarding the eigenstates 
and even analytic determination of the ground state in 
certain parameter regions |19|. Motivated by the pre- 
ceding, the main goal of this work is to provide an addi- 
tional careful study of the quadrumer spin-chain combin- 
ing complementary techniques and enhancing upon the 
parameter space investigated so far. 

The paper is organized as follows: In section |l|, we de- 
scribe the quadrumer-chain Hamiltonian more precisely 
and discuss the limit of decoupled quadrumers. Then we 
study t he i nter-quadrumer coupling Ji perturbatively in 
section III , We proceed with a numerical analysis of the 



model in section IV and try to fit the neutron scattering 
data with this method in section M. Finally, we summa- 
rize and discuss our results in section VI. Supplemen- 



tary results regarding the magnetization process of the 
quadrumer-chain model are presented in an appendix. 



antiferromagnetic, i.e. the physical regime is that with 
Ji, J2 > 0. In contrast, the physical regime for J' is less 
clear. Previous investigations of the model (Q) have as- 
sumed J' > ||l^,|l^ while, e.g., in [Q a ferromagnetic 
J' of approximately — 18meV has been used, focusing 
however on a different, i.e. zig-zag, charge order. We will 
therefore investigate both signs of J'. 

As a first step, we consider the case of decoupled 
quadrumers, Ji = where the Hamiltonian p) decom- 
poses into those of individual quadrumers (0). Rewriting 
the latter as 
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Tb,x = Six+2 + Six+i , 
Lx = Ta,x + Tb,x , 



(4) 



one infers the spectrum for an isolated quadrumer in Ta- 
ble |. The case J' > was discussed in detail in |l^ . The 
ground state is a singlet with a gap to S' = 1 excitations 
B for sufficiently large J' (J' > J2/2). For an illustration 
of the spectrum see [ [L5[ . The presence of a gap, however, 
is not immediately clear in two other regimes, namely 
for J' < — J2 where the ground state carries Lx = 2 and 
for — J2 < J' < J2/2 where it is degenerate between an 
Lx = and Lx = I representation. 



II. PRELIMINARIES 

The S = 1/2 Heisenberg Hamiltonian corresponding 
to the model in Fig. reads 
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A more conventional parameterization of the dimerized 
linear chain embedded in Fig. || would be Ji = (1 — 5)J, 
J2 = {1+S)J. This emphasizes its relation to the quarter- 
filled two-leg ladders of the high-temperature phase fin] 
whose low-energy properties are equivalent to spin-1/2 
Heisenberg chains. The dimerization is caused by the lat- 
tice distortions accompanying the phase transition [[7|,y. 
Consequently, estimates for the coupling constant of an 
effective Heisenberg chain such as J w 529K ?» 45meV 
p2[ can be used in the high-temperature phase. Such 
estimates are not expected to change drastically below 
the ordering temperature, apart from the dimerization. 
In particular the coupling along the chain should remain 
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TABLE I. Spectrum of a single quadrumer. 




FIG. 3. Effective model for the lowest magnetic excitations. 
Short arrows denote S = 1/2 spins and long ones S = 1. 
Coupling constants are indicated by the style of the lines. 

Even if the coupling Ji is switched on, Tb,x remains 
a good quantum number locally, i.e. for all x indepen- 
dently. This is related to the fact that local reflections 
of a single quadrumer at the middle chain are symme- 
try operations of the model in Fig. |[ Numerical re- 
sults (see sections IV and M below) show that the lowest 
magnetic excitations are always in the sector with all 



Tb,x = 1- These states are symmetric with respect to 
all local reflections of quadrumers along the chain. This 
finding yields important simplifications in determining 
the magnon-dispersion since it suffices to retain only the 
S = 1 combination of the two exterior S ^ 1/2 spins of 
a quadrumer. This leads to an equivalent description in 
terms of the simplified model sketched in Fig. ||. 



III. PERTURBATION THEORY IN Ji 

In this section we focus on a perturbative treatment 
of inter-quadrumer coupling. According to the previous 
section, there are three regimes with L^ = 0, 1 and 2 
ground states on each quadrumer. While the gap can be 
treated by standard non-degenerate perturbation theory 
in the case of L^ = 0, one has to resort to degenerate 
perturbation theory in all other cases. Here we show 
that this leads to effective Hamiltonians which also have 
a gap to magnetic excitations. 



A. J2/2 < J' 

For Ji — and J2/2 < J', the ground state carries 
Ta,x — Tb,x ~ ^, Lx = 0. In order to simplify the 
discussion, we will also treat J2 perturbatively, i.e. we 
will consider the regime J' ^ | Ji | , | J2 1 . Then the lowest 
excited state is obtained by creating an Ta,xo = Tb,xo = 
1, Lxg = 1 excitation on one quadrumer Xq. From this 
we obtain the following third-order expansion in Ji, J2 
of the dispersion 
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(5) 
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Eq. (|5|) agrees with the second-order result derived for 
J2 = in |lq|. One can read off from (ph that the min- 
imum of /Mk) is &i k — t: for ji > and at /c = for 
Ji < 0. 

Substituting a specific value for k into dq) , we can com- 
pute two further orders. At fc = and k = tt we then 
find the following fifth-order series for the gap: 
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The series presented above will be compared to numerical 
results in section K^ and one will see that they are quite 
accurate in the region where they are valid 123] . 



B. 



-J2 <J'< J2/2 



For — J2 < J' < J2/2, one reads off from Table | that 
there are two degenerate ground states on each cluster for 
Ji = 0. Both of them carry Ta,x = 0, one is a singlet with 
Tb,x = and the other one is a triplet with Tb.x = 1- 

In view of the effective spin-model of Fig. ^ only the 
triplet state with Tb,x = 1 contributes to the low-energy 
magnetic excitations for Ji / 0. Based on second-order 
degenerate perturbation theory in Ji one therefore ex- 
pects the system to be equivalent to an effective S* = 1 
chain as was already remarked in |15 1^. However, this 
effective chain is not a simple S — 1 Heisenberg-chain: 
While first-order interactions vanish, we find already at 
second order biquadratic terms which have not been dis- 
cussed before. Thus, the effective Hamiltonian is 
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x=\ *- 
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where Sx is an effective S* = 1 operator for the quadrumer 
X. We have determined the coupling constants in (0) to 
be 
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12 (J2 - 3 J') Jl + 17 (J' + Ji) J2J' - 8 J'* 
"32 (J2 - 2 J') (2J2 - J') (2J2 - 3 J') {.h + J') J2 ' 
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"2(J2- J')(4J2- J')(>/2-2J') ■ 



(10) 



These effective coupling constants are shown in Fig. for 
the region where the mapping applies. One finds that J 
is always positive (antiferromagnetic) while [3 is always 
negative. At the boundaries of the figure, we have the 
following divergences: J7 — > cxd as J' — > — J2 and (3 — > 
— cx) as J' -^ J2/2. These divergences are expected since 
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FIG. 4. Coupling constants for the effective biquadratic 
5=1 chain: The hnes denote J x J2 (full) and P (dotted). 



then leads to a small, but non-zero gap of the quadrumer 
model 

A « O.OO55J1 + 0{Jl) (12) 

for J' < — J2 and weakly antiferromagnetic Ji > 0. 

IV. LANCZOS DIAGONALIZATION 

We now proceed to study the spin gap over a wider 
range of parameters using Lanczos diagonalization. We 
have used periodic boundary conditions and mostly 
worked with the original model (|^), only for the case 
of L = 32 spin-1/2 spins have we resorted to the effective 
Hamiltonian of Fig. to reduce the dimensionality of the 
Hilbert space. 



they signal the limit of validity of the mapping to an 
effective S = 1 chain. 

The effective Hamiltonian (W with coupling constants 
( p^ ) is gapful with the exception of J' = and J' « 
O.4466J2 (corresponding io (3 — —1 - see for instance 
Sec. 6 of ^Jl). The ground state is the Haldane phase 
for — J2 < J' ^ O.4466J2 and it is spontaneously dimer- 
ized only in the small window 0.4466J2 ^ J' < >/2/2 
(see Sec. 6 of |Q). Interestingly though, the latter win- 
dow includes the region presumably most relevant to a'- 
NaV205, which we will find to be J' « J2/2 (see section 
M below). The biquadratic term thus has the important 
consequence that the relevant parameter region of the 
quadrumer chain is not adjacent to a Haldane state, but 
to a dimerized state. 

Regarding the size of the gap obtained by this map- 
ping, one observes that ^7 is very small in most of the 
region covered by Fig. El For example, for — 0.7264J2 < 
J' < J2/2, we have J2 x J^ < 1/40. In combination with 
the additional factor Jf this leads to an extremely small 
gap of the original Hamiltonian at least in the region 
Ji <C J2 where the perturbative approach is valid. 



C. J' < -J2 

Finally, for J' < — J2 and Ji = 0, the ground state 
of each quadrumer is an L^; = 2 representation. Us- 
ing first-order degenerate perturbation theory in Ji, the 
quadrumer chain can then be mapped to an effective 
S — 2 Heisenberg chain 
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(5=2) _ 
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L/4 



x — 1 



+1 
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where now Sx is an effective S = 2 operator for the 
quadrumer x. The gap of the effective S ^ 2 chain was 
estimated by DMRG to be Aoff. = (0.0876 ± 0.0013) Joff. 
[Hi. Substitution of Jeff. := Ji/16 into this DMRG result 
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FIG. 5. Spin gap as a function of J = Ji = J2 for J' > 0. 
Lines are for L = 8 (dotted), L = 16 (dashed) and L — 24 
(full). The symbols '-I-' denote an extrapolation to the ther- 
modynamic limit and the dashed-dotted lines [2,3] and [3,2] 
Fade approximants to the series (R) for J < and (H) for 
J>0. 

We begin with parameters that permit a comparison 
with perturbative results of the previous section. Fig. @ 
shows the spin gap, i.e. the minimum of the one-magnon 
dispersion, as a function of J for J' > and a homo- 
geneous chain with J = Ji = J2, including the region 
\J\ -^ J' where perturbation theory should be accurate. 
In this figure, we show numerical data for L = 8, 16 and 
24. Actually, we have determined the lowest excitations 
in the S"^ = 1 sector, but we are confident that this is 
in fact always an excitation with total spin S = 1. The 
finite-size data was extrapolated with a Shanks transfor- 
mation |p9|, yielding the pluses in Fig. H. These extrapo- 
lations are indistinguishable from the L = 24 data, show- 
ing that finite-size effects are extremely small at least for 
the parameters covered by the figure. 

Recall that eq. (0) shows that the minimum of the one- 
magnon dispersion is located at fc = for J < and at 
k = -K for J > 0. This is confirmed by our numerical data. 



Therefore, one should compare the numerical data to the 
series (0) and (^) for J < and J > 0, respectively. 
Pade approximants to these series are shown in Fig. 
and one observes good agreement for \J\ < J'. While 
the approximants stay close to the numerical data for 
negative J as large as J = —2.5 J', systematic deviations 
are observed at large positive J. This can be attributed 
to the fact that other S = 1 excitations start to mix in 
at J' w J (compare Table | and [|l6|). 

As a first summary, we find good agreement between 
our numerical data and perturbation theory in the regime 
where the latter should be accurate. For J ^ J', however, 
the numerical approach is far superior since the nature of 
the lowest 5 = 1 excitation becomes more complicated. 
Yet, finite-size effects still remain small. 
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FIG. 6. Same as Fig. Bl but now as a function of J' for 
J > 0. The symbols 'x' denote L — 32 data and a ' + ' with 
error bars an extrapolation to the thermodynamic limit. 



gion — 2J ^ J' ^ —J. In fact, while A/J is already of 
the correct order of magnitude, it still has to decrease 
substantially for J' — > — cx) in order to reach the value 
given by (|l^). In the limit J' -^ — oo, we also have to 
recover the correlation length ^ of the 5* = 2 chain which 
is known to be ^ w 50 clusters {i.e. L w 200 - see for 
example p7|), thus leading to the increasing finite-size 
effects in Fig. g as J' goes to larger negative values. 

In the region of positive J, J' > 0, our results can be 
compared to DMRG resuhs (see Fig. 2 of ||l|l). While 
we find agreement for large J', deviations of the order 
of J/50 can be observed in the region of small J'. For 
instance, we extrapolate a gap A « 0.077J for J' — 0.5 J 
while the corresponding value o f llSll is A w 0.05 J. Fur- 
thermore, the DMRG result of ||15[| for the gap is essen- 
tially zero already for J' = 0.4J while our extrapolation 
at J' = 0.35J yields A = 0.0201(6)J which is clearly 
non-zero. This discrepancy can be observed already at a 
fixed system size (e.g. L — 32) [2^], i.e. it is not due to 
the extrapolation. There are two main distinctions be- 
tween our Lanczos and the DMRG results of jl^]: First, 
DMRG is subject to truncation errors which are absent in 
our approach. Second, we use periodic boundary condi- 
tions in contrast to open ones in p5| . Periodic boundary 
conditions typically have smaller finite-size effects than 
open ones and open boundary conditions can even lead to 
further boundary excitations. We therefore believe that 
extrapolation of our Lanczos data for system sizes up to 
L = 32 yields the most accurate results obtained so far 
for the thermodynamic limit of the model (j^ . 



V. FITTING THE NEUTRON SCATTERING 
DATA 



We now concentrate on the region J > 0, i.e. the one 
which is appropriate for ck'-NaV205, and study the de- 
pendence on J', focusing on small and negative values. 
Fig. 1^ again shows the spin gap, now as a function of J' 
for J = Ji = J2 > 0. The extrapolations have again 
been performed with a Shanks transformation [Q. The 
actual estimate is based on the L = 16, 24 and 32 data 
with error bars determined from the difference with re- 
spect to an extrapolation based on i = 8, 16 and 24. For 
J' > 0, the trend observed in Fig. |[ i.e. that L = 8, 16 
and 24 is completely sufficient for the extrapolation, is 
confirmed. For J' < and small, the errors first remain 
small with an increasing uncertainty as J' becomes large 
and negative. Larger error bars indicate larger finite-size 
effects which in turn originate from an increasing corre- 
lation length. 

Fig. 1^ shows that the gap vanishes at J' = at any fi- 
nite system size, the reason being that the exterior spins 
are free in this case and give rise to excitations with ex- 
actly zero energy. For any non-zero J', in contrast, the 
gap seems to remain non-zero even in the thermodynamic 
limit. 

A maximum of A/J emerges in the ferromagnetic re- 



We have seen so far that there are several regions where 
different mechanisms open a spin gap in the quadrumer- 
chain model with a large variation of its size. While this 
looks promising for the applicability of the model to a'- 
NaV205, the real test is whether the neutron-scattering 
results for the dispersion along the 6-axis |13,M can be 
fitted with reasonable parameters. 

Here we follow |13] and use the foUowingcriteria: (i) 
The low-temperature experimental data |12| , |l4| has a 
miminum in the center of the high-temperature Brillouin 
zone. Taking into account doubling of the unit cell (as 
is appropriate for the quadrumer-chain model), the min- 
imum of the dispersion should therefore be located at 
fc = 27r = and should have a value A(0) « lOmeV. (ii) 
Since the magnon dispersion was recentely traced up to 
energies of 40meV jlj], one should have A(7r) > 40mcV. 

In Figs. Ra),b) we first show the dispersion of the low- 
est S — I excitation corresponding to the two points 
J' = ± J/2 in Fig. 0. One observes that the minimum is 
located at fc = tt and remains there after extrapolation 
to the thermodynamic limit. Fig. ^) shows the disper- 
sion for J' = 0.65 J2 and a small dimerization S « 0.053 
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FIG. 7. Dispersion of lowest 5* = 1 excitation for a) 
Ji = J2 = 2J', b) Ji = Ja = -2 J' and c) Ji = 0.9 J2, 
J' — O.65J2. Symbols are for i = 16 (open circles), L = 24 
(filled circles) and L = 32 (pluses). Lines are guides to the 
eye. 



which is close to a parameter set proposed as a possible 
fit for a'-NaV205 in ||l^. One observes that our data is 
already well converged with system size and one can thus 
easily extrapolate it to L = 00. We find A(0) ~ 0.239 J2, 
which agrees with [|l6| , but the minimum is still at /c = tt: 
A(7r) « O.215J2, leading to A(7r)/A(0) « 0.9. We there- 
fore agree with [|5| that this parameter region is not ap- 
propriate for a'-NaV205. 



A(0)/A(7i) 
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FIG. 8. A(0)/A(7r) as a function of coupling to exter- 
nal spins J' / J2 and coupling between quadrumers J1/J2 for 
L = 16 (dotted lines) and L = 24 (full lines). 

In order to test whether a good fit can be obtained 
at all, we have performed a systematic numerical scan 
of A(0)/A(7r) at L = 16 and L = 24. The range of 
interest including ferromagnetic couplings J' < we find, 
is depicted in Fig. 0. Apart from this range we have also 
studied a wide region of J'/ J2 and Ji/ J2. We found little 
variation outside the parameter window shown in Fig. 0. 
Note that A(fc) was determined for S^ = \ excitations 



which do not necessarily carry 5=1, but may have 
higher spin. Furthermore, there should be a two-particle 
state at fc = with an energy not larger than 2A(7r) in 
the thermodynamic limit. The ratio A(0)/A(7r) should 
therefore not exceed a value of 2 in the thermodynamic 
limit unless the gap closes. This indicates large finite 
size effects in some regions of Fig. @. However, we are 
interested in parameters where A(0)/A(7r) < 1. In that 
case no complications due to multi-magnon excitations 
are expected. 
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FIG. 9. A(0)/A(7r) for Ji = J2 = J as a function of J' 
for 1/ = 16 (open circles), L = 24 (filled circles) and L = 32 
(pluses). Lines are guides to the eye. 

For clarity. Fig. S shows a section along Ji = J^ of 
Fig. bI including some L = 32 data points (compare also 
Figs.1^),b)). One observes that finite-size effects can be 
important, in particular for J' < where A(0)/A(7r) 
exceeds the aforementioned limiting value of two. It 
is evident from this figure that finite-size effects reduce 
A(0)/A(7r). The condition A(0) < A(7r) appears nev- 
ertheless impossible to satisfy at J\ — J^- Moreover, in 
the entire region with J' < 0, Fig. g clearly shows that 
A(0) is always larger than A(7r), restricting the possible 
parameters for a fit to lie in the region J' > 0. This is 
consistent with refs. |13,M. 

On the other hand, there is a narrow region around 
J' = J2/2 with J\ < J2 where A(0) becomes smaller 
than A(7r). An example of an excitation spectrum in 
this region is shown in Fig. O. One does indeed see 
a single magnon excitation with a minimum at fc = 
whose dispersion is interpolated by the full line. This 
dispersion would be consistent with early neutron scat- 
tering data p^ . However, very recent measurements |l4] 
clearly show that the bandwith in Fig. |l0| is too small. 
One may obtain a larger bandwith within the present 
model by changing the coupling constants slightly, how- 
ever, only at the expense of reducing A(0). In order to 
have A(0) = lOmeV and A(7r) > 40meV one then needs 
a coupling constant J — J2/(l + 6) }^ 300meV with a 
dimerization^ w0.38...0.54and J' « (0.5 .. .0.55) J2 [|9). 
Both, the large values of J and of the dimerization are 
certainly too large to be plausible (see section O) . There- 
fore, even though our numerical results differ in detail we 
agree with 1l^ regarding that the cluster spin model dll) 
does not result in a quantitive fit for a'-NaV205 with 
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FIG. 10. Low-lying excitations at J1/J2 ~ 7/20, 

J' /J2 = 11/20. Open symbols denote data for L — 16 and 
filled symbols for L — 24. Boxes stand for S — excita- 
tions, circles for S — 1 excitations and triangles for S — 2 
excitations. The pluses show the lowest S — 1 excitations 
for L = 32 at fc = and tt. The full line denotes the lowest 
S — 1 excitation (magnon) whereas the dashed line shows the 
lowest, non-dispersive S — excitation. The shaded area in- 
dicates the expected location of the two-magnon continuum. 



plausible parameters. 

Among the remaining excitations displayed in Fig. |lC|, 
there is a localized (fc-independet) singlet with A(fc) « 
O.I4J2 (indicated by the dashed line). This excitation can 
be interpreted as a singlet Tb.xo = in one quadrumer 
xq (while all other Tb,x — 1) which is prevented from 
propagating since Tb,x is locally conserved. Note that 
such excitations are not present in the effective model of 
Fig. ||, but only in the original one of (|l|). 

We would like to conclude this section with a remark 
on Fig. H: There is one dip in A(0)/A(7r) for Ji < J2 both 
for negative and positive J'. These dips become sharp 
for Ji — > where they signal transitions between different 
local ground states. In this limit, the dips are located at 
J' /J2 = —1 and J'/J2 = 1/2 where according to Table || 
the transitions occur at Ji = 0. In general, our numerical 
data is consistent with a vanishing gap at the location of 
these dips, i.e. a continuous transition between different 
phases. However, due to the large correlation lengths 
expected in some of the phases, larger system sizes would 
be needed in order to determine the phase diagram of the 
model (|l]) accurately. 



VI. CONCLUSIONS AND OUTLOOK 

We have studied an S* = 1/2 quadrumer-chain model 
using perturbation theory and the Lanczos method. This 
model has been proposed to explain the magnetic excita- 
tions of a'-NaV205 in its low-temperature phase |^,g,|l5[. 
We have shown the quadrumer-chain Hamiltonian to dis- 



play a rich phase diagram which we believe the present 
study has just begun to unveil (compare Fig. g): Differ- 
ent phases are smoothly connected to the local Lx — 0, 1 
and 2 ground states at Ji = 0. In addit ion, t he mapping 
to an effective S — 1 chain in section IIIB shows that 
the Lx = 1 region consists of a Haldane and a dimer- 
ized phase at least at small Ji. The latter phase arises 
because of a biquadratic interaction which has not been 
realized in previous studies. However, further work is 
needed to determine all of the phase boundaries of the 
quadrumer chain accurately. 

The model (|l|) gives also rise to interesting behavior in 
an external magnetic field. The case Ji = J2 is discussed 
in an appendix where we show that in addition to a spin 
gap one also finds a plateau in the magnetization curve 
at half the saturation ma gne tization. This is similar as 
for the model studied in |21||, though it remains to be 
investigated whether one can also obtain other plateaux 
in the present model for Ji ^ J2 like the one at a quarter 
of the saturation magnetization found in 12^ . 

Finally, we have assessed the relevance of the 
quadrumer-chain model to the magnetic excitations of 
a'-NaV2 05. We found that in order to fit the neutron 
scattering data for the dispersion along the chain di- 
rection jlj] one has to resort to parameters which are 
not very plausible. A good fit can even be obtained for 
Ji < 0, (compare (||) and Fig. ||) although ferromag- 
netic exchange along the chain is clearly an unphysical 
choice. The applicability of the quadrumer model to a'- 
NaV2 05 may therefore be questioned [|l5[, leaving the 
proper microscopic model for this material an open is- 
sue. In fact, other models proposed in this context are 
similarly deficient. Assuming, e.g., pure zig-zag charge 
order [g[, the lowest S = 1 excitation would belong to 
a two-particle continuum (see, e.g., |30|) rather than the 
experimentally observed magnon state |lj,|lj]. A possi- 
ble remedy is to attribute the spin gap to dimerization 
pll rather than frustration. However, this proposal has 
been disputed |1^] based on the spatial symmetries of the 
charge ordered state. Clearly, further experimental input 
is needed to decide these issues. For the moment, we can- 
not rule out that the correct description of a'-NaV205 
will turn out to be a modification of the cluster model 
discussed in the present work. 
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APPENDIX: MAGNETIZATION PROCESS 

In this appendix we discuss the magnetization process 
of the quadrumer-chain model at Ji ~ J2 ~ J- AppU- 
cation of an external magnetic field amounts to adding a 
term — /iX]£c=i ^x ^^ the Hamiltonian (|l|). We have com- 
puted the magnetization curve numerically for i < 32 
(as one example, the magnetization curve for J — J' is 
shown in the inset of Fig. Hi]) and then extrapolated it 
to the thermodynamic limit. The final result is shown 
in Fig. O. First, we find an (Af) = plateau which 
is equivalent to the spin gap (accordingly, its boundary 
curve in Fig. O is equivalent to Fig. ra). In addition, we 
find a clear plateau at (A/) = 1/2 (m a normalization 
where the fully polarized state has (M) — 1). Such a 
plateau is expected on the basis of the limit of decoupled 
quadrumers (g) - see ||32|] and references therein. Simi- 
lar behavior has been observed for the related model of 
ref. [Il9| where an interesting sequence of magnetization 
plateaux was found [pl|. One remarkable feature of the 
present model is that there is evidence for a transition in 
the (M) = 1/2 plateau state at J « 0.65 J'. 




FIG. 11. Groundstate phase diagram in a magnetic field 
h for Ji = J2 = J. White areas denote parameter regions 
with magnetization plateaux with the indicated values of the 
magnetization (Af ) whereas the grey shaded regions denote 
smooth transitions. Inset: Magnetization curve at J' = J for 
L — 32 (thin line) and an extrapolation to the thermodynamic 
limit (bold line). 
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